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Abstract 

Promptly decaying lightest neutralinos and long-lived staus are searched for in 
the context of light gravitino scenarios. It is assumed that the stau is the next 
to lightest supersymmetric particle (NLSP) and that the lightest neutralino 
is the next to NLSP (NNLSP). Data collected with the Delphi detector at 
centre-of-mass energies from 161 to 183 GeV are analysed. No evidence of the 
production of these particles is found. Hence, lower mass limits for both kinds 
of particles are set at 95% C.L.. The mass of gaugino-like neutralinos is found to 
be greater than 71.5 GeV/c 2 . In the search for long-lived stau, masses less than 
70.0 to 77.5 GeV/c 2 are excluded for gravitino masses from 10 to 150 eV/c 2 . 
Combining this search with the searches for stable heavy leptons and Minimal 
Supersymmetric Standard Model staus a lower limit of 68.5 GeV/c 2 may be set 
for the stau mass independent of the mass of the gravitino. 
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1 Introduction 

In models including supersymmetry (SUSY), it is often assumed that the messengers 
of supersymmetry breaking couple to the observable sector with interactions of gravita- 
tional strength and that the SUSY breaking scale in the hidden sector is of the order 
of 10 11 GeV. An alternative possibility is that supersymmetry is broken at some lower 
scale (below 10 7 GeV), and that the ordinary gauge interaction acts as the messenger 
of supersymmetry breaking |l|J^]. In this case, the gravitino, G, is naturally the lightest 
supersymmetric particle (LSP) and the lightest Standard Model superpartner is the next 
to lightest supersymmetric particle (NLSP). Thus, the NLSP is unstable and decays to 
its Standard Model (SM) partner and a gravitino. 

Since the gravitino couplings are in general, with the exception of the so-called ultra- 
light gravitino scenarios, suppressed compared to electroweak and strong interactions, 
decays to the gravitino are in general only relevant for the NLSP and therefore the 
production and decay of supersymmetric particles at high energy colliders would genera- 
lly take place through Standard Model couplings Q The supersymmetric particles decay 
into the NLSP, which eventually decays to its SM partner and a gravitino. The specific 
signatures of such decays depend crucially on the quantum numbers and composition of 
the NLSP. 

Although most of the attention has been focused on the case where the neutralino is 
the NLSP, it is also possible that the NLSP is any other sparticle, and in particular a 
charged slepton. The number of generations of supersymmetry breaking messengers in 
minimal models, n, determines over most of the parameter space which particle is the 
NLSP [0,f||)]||. For example, for one generation of messengers, the lightest neutralino 
tends to be the NLSP, while for two or more generations, right handed sleptons are 
favoured. Moreover, when left-right sfermion mixing [|7| occurs, the corresponding f 
state, fi, becomes the NLSP. 

Throughout this work, it is assumed that the fi is the NLSP and that the lightest 
neutralino, Xii ls the next-to-NLSP (NNLSP). The f x width is given (independently of 
the f mixing) by the two-body equation: 

~ TO- 

m^T + G)= ;i 2 (i) 

where is the mass of the fi, vn^ is the mass of the G and M p is the Planck mass 
(2.4 x 10 18 GeV). In the last equation, the mass of the r has been neglected. The mean 
decay length obtained from equation ([]]): 

-5 / \ 2 

I m G 

.100GeV/cV UeV/c 

depends strongly on mf, ttlq and the energy of the f\, E. The dependence of the mean 
decay length, L, on ttlq could be also interpreted in terms of the supersymmetry breaking 
scale, yF, through the relation: 

mr= !- -2.5 | | eV/c 2 (3) 

° V3M P \100 TeV J ' K J 

For y/F < 1000 TeV (m^ < 250 eV), the decay can take place within the detector. This 
range of \fF is in fact consistent with astrophysical and cosmological considerations ||. 

x One exception to this rule being the process e + e — — > Z* /-/* —* Gxi f° r the case of ultra-light G scenarios. 



L = 1.76 x lO-^Vm^ - 1)2 n ^ ) ( ) cm ' ( 2 ) 
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Two searches are presented here. Firstly, Xi P a i r production with Xi decaying to f\T 
and f\ then decaying promptly into tG. The signature of the signal would be four r's 
with missing energy and momentum from the two gravitinos (in addition to the energy 
and momentum carried away by the neutrinos of the decay of the r's). 

The second search concerns f\ pair production followed by the decays fi — > tG within 
the detector volume. The signature of such an event will be a track of a charged particle 
with a kink or a decay vertex when the f\ decays inside the tracking devices. If the 
decay length is too short (small tUq) to allow for the reconstruction of the f% track, only 
the decay products of the r will be seen in the detector, and the search will then be 
based on track impact parameter. However, if the decay takes place outside the tracking 
devices (large mg), the signature will be that of a heavy charged particle already studied 
in DELPHI ||. For very light uiq the limits from the search for MSSM stau can be 
applied [fL0|j . All these searches have been combined to obtain a limit on vrif R independent 
of the G mass. 

The data samples and event selections are respectively described in sections |2| and 3, 
while the results are presented in section |j. It will be seen in section ^ that these two 
searches, together with those for Xi ~~ * lG \\ 1| (in the x? NLSP scenario) and promptly 



decaying f\ pair production JlOj complement each other for different domains of the 
gravitino mass. 



2 Event sample and experimental procedure 



The search for neutralino pair production is based on data collected by the DELPHI 
experiment during 1996 and 1997 at centre-of-mass energies of 161, 172 and 183 GeV. 
The total integrated luminosities for the three centre-of-mass energies are 9.7, 10.4 and 
53.9 pb _1 respectively. The search for stau pair production with big impact parameters 
and secondary vertices is based on data collected by the DELPHI experiment during 
1997 since the results obtained with the data collected in 1995 (at ^fs = 130 — 136 GeV) 
and 1996 are already published in ||12| . The present analysis updates those results. The 
search for stau pair production with small impact parameters is based on data collected 
from 1995 to 1997. A detailed description of the DELPHI detector can be found in fl3| 
and its performance in 
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To evaluate the signal efficiencies and background contaminations, events were genera- 



ted using different programs, all relying on JETSET 7.4 [15], tuned to LEP 1 data [|H| for 
quark fragmentation. The program SUSYGEN |7l was used to generate the neutralino pair 
events and their subsequent decay products. In order to compute detection efficiencies, a 
total of 3000, 10000 and 14000 events were generated with centre-of-mass energies of 161, 
172 and 183 GeV respectively, and masses 47 GeV/c 2 < m fl + 2 GeV/c 2 < m^o < y/s/2. 
A stau pair sample of 18000 events (subdivided in 15 samples) was produced with PYTHIA 
5.7 [0 at 183 GeV centre-of-mass energy, the staus having mean decay lengths from 
0.25 to 1000 cm and masses from 40 to 90 GeV/c 2 . Another sample of 35000 stau pairs 
was produced with SUSYGEN for the small impact parameter search (see below), with 
centre-of-mass energies ranging from 130 GeV up to 183 GeV. 

The background process e + e~ — > qq{n / ~f) was generated with PYTHIA 5 . 7, while DYMU3 
l"8| and K0RALZ JTIJ were used for fi + fi~(j) and t + t - (7), respectively. The generator of 
reference [2D] was used for e + e~ — > e + e~ events. 

Processes leading to four-fermion final states, (Z /'-/)* (Z / 'j)* , where * means of-the- 
mass-shell, W + W~, Weu e and Ze + e~, were also generated using PYTHIA. The calculation 
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of the four-fermion background was verified using the program EXCALIBUR pl |, which 
consistently takes into account all amplitudes leading to a given four-fermion final state. 

Two-photon interactions leading to hadronic final states were generated using TWOGAM 
22|| , separating the VDM, QPM and QCD components. The generators of Berends, 



Daverveldt and Kleiss [23] were used for the leptonic final states. 

The cosmic radiation background was studied using the data collected before the 
beginning of the 1997 LEP run. 

The generated signal and background events were passed through the detailed simu- 
lation |14] of the DELPHI detector and then processed with the same reconstruction and 
analysis programs used for real data. 



3 Data selection 

3.1 Neutralino pair production 

In this section, the selections used to search for the process e + e~ — > XiXi ~^ tittxt — ► 
tGttGt are presented. 

The reconstructed tracks of charged particles were required to have momenta above 
100 MeV/c and impact parameters below 4 cm in the transverse plane and below 10 cm in 
the longitudinal direction. The relative error on the measurement of the momentum was 
to be smaller than 100%. Clusters in the calorimeters were interpreted as neutral particles 
if they were not associated to charged particles and if their energy exceeded 100 MeV. 
All charged and neutral particles that satisfy these criteria are considered good particles 
and they are used to compute the relevant event quantities. To assure good quality 
of the data, the ratio of good to total number of tracks was required to be above 0.7. 
Tracks that did not pass quality selection but had an associated calorimetric energy of 
at least 2 GeV had their angles taken from those of the track, but their momentum was 
recomputed as that of the calorimetric measurement. Such tracks were not included in 
the good sample. Events had to have between four and ten good charged particle tracks. 
In addition, it was required that the thrust be less than 0.99. The transverse momentum, 
computed as the transverse component with respect to the beam axis of the vector sum of 
the momenta of good charged and neutral particles, pt, had to be bigger than 3 GeV/c. 
And the absolute value of the cosine of the polar angle of the missing momentum vector be 
less than 0.95. Very forward-going events were eliminated by requiring that the energy 
in a cone of 30°, E30, around the beam-pipe to be less than 70% of the total visible 
energy, E vis . With this preselection, the total number of simulated background events 
and real data events was reduced by a factor of about 6000. Only events passing these 
pre-selections were analysed further. 

The selection takes advantage of the fact that signal events can be separated into 
two different kinematic regions of the (77^0,7?%^) space: when the mass difference Am = 
m^o —m?! is bigger than about 10 GeV/c 2 , all four r's carry similar momenta. When the 
difference becomes smaller, the two r's coming from the decay of the fi tend to be the 
most energetic, increasingly so as the x? mass increases. The Durham algorithm [25] was 
used to divide the event in four jets by allowing y cut to vary as a free variable. Numbering 
the jets from 1 to 4 with E x > E 2 > E 3 > E 4 , a variable r was defined as: 

E3 + E4 

r = — — . (4) 

Ei + E 2 
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An example of the distribution of r for simulated samples with two values of Am can be 
seen in figure |l| It should be noticed that the distribution of r shifts towards lower values 
with increasing neutralino masses. The simulated background samples were then divided 
into two samples above and below r = 0.1 and different requirements were imposed in 
the two cases. 
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Figure 1: Two examples of the distribution of the variable r (see text). The positive-slope 
hatched histogram shows r for Am = 3 GeV/c 2 . The negative-slope hatched histogram 
shows r for Am = 20 GeV/c 2 . 



Two sets of cuts were applied in order to reduce the 77 and ff (7) backgrounds and a 
third set of cuts to select events according to their topology: 

1- Cuts against 77 backgrounds: the transverse energy, Et, should be bigger than 
11 GeV for r > 0.1 (Et > 12 GeV for r < 0.1). The energy in a cone of 30° 
around the beam axis was further restricted to be less than 60% of the total visi- 
ble energy to avoid possible bias from the Monte Carlo samples. The missing mass 
should be smaller than 0.88-y/s (0.9a/s). The momentum of the charged particle with 
largest momentum should be bigger than 4 GeV/c (3 GeV/c). These cuts reduced 
the 77 background by a factor of the order of 30. 

2- Cuts against ff(7) backgrounds: the number of good tracks should be smaller than 
7 (9). The maximum thrust was further reduced from 0.99 to 0.975. Dividing each 
event into two jets with the Durham algorithm, its acoplanarity should be bigger 
than 8°. The missing mass of the events should be bigger than 0.3^. After these 
cuts, the ff (7) background was reduced by a factor of the order of 15. 

3- Cuts based on topology: signal events tend naturally to cluster into a 4-jet topology. 
All jets should be at least 17° away from the beam direction. When reduced by the 
jet algorithm into a 2-jet configuration, the charged particles belonging to each of 
these jets should be in a cone broader than 20°. Finally, the axes of each of the four 
jets should be separated from the others at least by 8° (4°). 
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Figures |2| to [| show some of the distributions relevant for these selection criteria at 
y/s = 183 GeV. Table [1] shows the effect of these cuts at y/s = 183 GeV on the data, 
expected background and the signal for m^o = 75 GeV/c 2 and m fl = 55 GeV/c 2 . The 
discrepancy between data and simulation on the last bin of figure |2|-a is attributed to the 
poor description of 77 events in the simulation. 



Cut 


77 


ffl 


4-fermion 


Total MC 


Data 


Signal 


pre-selection 


496 ± 16 


44.5 ± 1.5 


13.1 ±0.6 


554 ± 16 


567 


61.4% 


1 


18 ±2 


40.6 ± 1.4 


12.1 ±0.6 


70.8 ±2.6 


84 


59.2% 


2 


2.2 ±0.6 


2.9 ±0.4 


4.5 ±0.4 


9.6 ±0.8 


12 


45.4% 


3 





0.23 ±0.09 


0.27 ±0.07 


0.50 ±0.11 


2 


38.3% 



Table 1: Number of events remaining in the data and simulated samples at y/s = 183 GeV 
after various stages of the selection procedure described in the text. The signal efficiencies 
corresponds to m^o = 75 GeV/c 2 and m fl = 55 GeV/c 2 . 

After these cuts, an efficiency between 25 and 45% was obtained for the signal events, 
and estimated backgrounds of 0.12 ± 0.08, 0.15 ± 0.09 and 0.50 ± 0.11 events at y/s = 
161, 172 and 183 GeV, respectively. 



3.2 Stau pair production 

This section describes the selection criteria used in the search for the process e + e~ — > 
fifi — > tGtG. As described in section |], the mean life-time of the f\ depends on the 
mass of the gravitino. Thus, for a gravitino with a mass of the order of a few hundred 
eV/c 2 or more, the stau would be sufficiently long lived to decay outside the detector. 
When the mass of the gravitino is between a few eV/c 2 and a few hundred eV/c 2 , one 
or both staus would decay in flight in some part of the detector, creating a well defined 
secondary vertex. The search for these decays is described in subsection |3.2.1| . If the mass 
of the gravitino is even smaller, stau pair production would produce displaced vertices. 
This search is described in subsections |3.2.2| and |3.2.3| . 

3.2.1 Search for secondary vertices 



This analysis exploits a peculiarity of the f\ — > tG topology in the case of intermediate 
mass gravitinos, namely, one or two tracks coming from the interaction point and at least 
one of them with either a secondary vertex or a kink. 



Rather loose preselection cuts, similar to those presented in reference [|12| were imposed 
on the events in order to suppress the low energy background (beam-gas, beam- wall, etc), 
77, e + e~ and hadronic events. To compute the following quantities the reconstructed 
tracks of charged particles were required to have momenta above 100 MeV/c and impact 
parameters below 4 cm in the transverse plane and below 10 cm in the longitudinal 
direction. Clusters in the calorimeters were interpreted as neutral particles if they were 
not associated to charged particles and if their energy exceeded 100 MeV. However, no 
quality requirements were imposed on the reconstructed tracks in the following stages. 

• Charged particle multiplicity between 1 and 10; 

• visible energy above 10 GeV; 
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P.eadi„ g (GeV/c) 



Figure 2: (a) missing mass and (b) momentum of the leading charged particle, for 
data (dots), Standard Model simulation (cross-hatched histogram) and one of the sim- 
ulated signals with cross-section not to scale (blank histogram) after preselection at 
\/s = 183 GeV. The arrows indicate selection criteria imposed as explained in the text. 
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Figure 3: Acoplanarity of data (dots), Standard Model background simulation 
(cross-hatched histogram) and one of the simulated signals with cross-section not to 
scale (blank histogram) at y/s = 183 GeV, after the cut to remove 77 events. The arrow 
indicates selection criterion imposed as explained in the text. 



• total electromagnetic energy below 40 GeV; 

• transverse momentum, computed as the transverse component with respect to the 
beam axis of the vector sum of the momenta of charged and neutral particles, pt, 
greater than 5 GeV/c; 

• energy measured in the very forward calorimeters below 10 GeV. 

These preselection cuts leave about 0.6% of the whole data sample. 

The tracks of the events that survived the preselection cuts were grouped in clusters 
according to their first measured point (starting point). This clustering procedure is 



described in |[26|| . Each cluster contained all tracks whose starting points differ by less 
than 2 cm. The starting point of a cluster was defined as the average of the starting points 
of its tracks. This procedure allowed for clusters with a single track if its momentum was 
larger than 1.5 GeV/c. Events were rejected if more than 6 tracks were not grouped 
in clusters or no cluster was found. A cluster with only one track was considered a fi 
candidate track if its trajectory was compatible with that of a particle coming from the 
interaction point (according to the selection criteria described in reference [^]) and its 



momentum was greater than 2 GeV/c. 

For each such fi candidate, a search was made for a second cluster with starting point 
radius in the transverse plane (xy plane) greater than that of the first measured point of 
the track of the fi candidate, and an angular separation between the directions defined 
by the beam spot and the cluster starting points smaller than 90° in the xy plane. This 
second cluster was assumed to be formed by the decay products of the r coming from the 
fi — ► tG process. Therefore, the f± candidate and the r cluster had to define a secondary 
vertex. If the r cluster included more than one charged particle, only the one with the 
highest momentum was used to search for the decay vertex or kink (crossing point with 
the fx track). 
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Figure 4: (a) minimum angle of a jet to the beam, (b) maximum angular broadness 
of a jet at the 2-jet level and (c) angle between jets, for data (dots), Standard Model 
background simulation (cross-hatched histogram) and one of the simulated signals with 
cross-section not to scale (blank histogram) at \/s — 183 GeV, after the cut to remove 
ffil) events. The arrows indicate selection criteria imposed as explained in the text. 
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The tracks were parametrised with respect to their perigee |27j to calculate the point 
of closest approach between the two tracks (the candidate f\ track and the selected track 
from the candidate r cluster). The conditions to define a good crossing point between 
the track of the f\ and the selected track of the r decay candidates are described in 



reference [12] 



Fake decay vertices could be present amongst the reconstructed secondary vertices, 
being produced by particles interacting in the detector material or by radiated photons 
if the particle trajectory was reconstructed into two separated tracks. To eliminate these 
classes of events, additional requirements were imposed: 

• to reject hadronic interactions, any reconstructed hadronic interaction (secondary 
vertices reconstructed in region where there is material) must be outside a cone of 
half angle 5° around the kink direction; 

• to reject photon radiation in the case of r clusters with only one track, there had 
to be no neutral particle in a 3° cone around the direction defined by the difference 
between the f\ momentum and the momentum of the r daughter calculated at the 
crossing point; 

• to reject segmented tracks, the angle between the tracks used to define a vertex had 
to be larger than 6°. 

If no pair of tracks was found to survive these conditions, the event was rejected. Fi- 
gure [| shows the distribution of these three angles for real data, expected Standard 
Model background simulation and simulated signal for = 60 GeV/c 2 decaying with 
a mean distance of 50 cm. The excess of data in the first bins of figure [5] (c) is due to 
underestimation in the simulation of mismatchings between the tracking devices. 

One event in real data was found to satisfy all the conditions described above. The 
event was the superposition of a low energy event with a cosmic muon crossing the 
detector. However, the two tracks of the cosmic muon follow the cosmic muon rejection 
criteria used on subsection [3.2.2| based on impact parameter and on acollinearity. Thus, 



the event will not be considered as a candidate. This kind of events was not simulated 
and therefore its removal does not affect the calculated efficiencies. 

The vertex reconstruction was sensitive to decay lengths in the xy plane, R, between 
15 cm and 90 cm. Within this region a vertex was reconstructed with an efficiency 
of ~54% since the VD (Vertex Detector) and the ID (Inner Detector) were needed to 
reconstruct the f\ track and the TPC (Time Projection Chamber) to reconstruct the 
decay products. The efficiency is flat inside the sensitive region and drops to zero for f\& 
decaying near the outer surface of the TPC. The shape of the efficiency distribution was 
independent of the f\ mass; it simply scaled down near the kinematic limit. The loss of 
efficiency near the kinematic limit is due to the fact that the f\ boost is smaller and the 
vertex reconstruction less efficient when the angles between the f\ and the r products 
increase. 

The efficiencies for different mean decay lengths and f\ masses were calculated by 
applying the above selections to the simulated signal samples. For a 60 GeV/c 2 f\ with 
mean decay length of 50 cm the vertex search efficiency is of the order of 55%. 

3.2.2 Large impact parameter search 



To investigate the region of low gravitino masses (short decay lengths) the previous 
search was extended to the case of the f\ decaying between 0.25 cm and around 10 cm. 
In this case the fi track was not reconstructed in the ID and only the r decay products 
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Figure 5: (a) angle between the hadronic interaction and the reconstructed vertex, (b) 
angle between the electromagnetic shower and the direction defined by the difference 
between the momenta of f\ and its associated r, defined at the crossing point, and (c) 
angle between the tracks of the kink, for real data (dots), expected Standard Model back- 
ground (cross-hatched histogram) and simulated signal for = 60 GeV/c 2 decaying 
with a mean distance of 50 cm (blank histogram). The arrows indicate selection criteria 
imposed as explained in the text. 
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were detected. The events used in this search contained exactly two single track clus- 
ters (i.e. two charged particles with momentum larger than 1.5 GeV/c and a distance 
between starting points greater than 2 cm) which were acollinear and had large impact 
parameters^]. Cosmic rays, badly reconstructed tracks or interactions in the detector ma- 
terial could result in large impact parameters. However, the two tracks in a cosmic event 
usually have impact parameters of the same order and opposite sign. The acollinearity 
in back-to-back events with badly reconstructed tracks or interactions was always small. 
Figure |5] shows the scatter plot of the maximum impact parameter versus the minimum 
one in the R<f> plane. Figure ^ shows the acollinearity distribution for events with two 
tracks in the TPC. Simulated signal events with = 60 GeV/c 2 and a mean decay 
length of 2.5 cm are compared with cosmic muon events, simulated Standard Model back- 
ground and real data. The data points in figure [7] contain cosmic radiation events that 
are not simulated. 




Figure 6: Impact parameters of two track events for a simulated signal of 
= 60 GeV/c 2 with mean decay length of 2.5 cm, cosmic muons, Standard Model 
expected background and real data. It is plotted the maximum impact parameter versus 
the minimum one in the R(j) plane. The area between the lines was excluded by the 
cosmic rejection criteria as described in the text. 

The impact parameter search was only applied to those events accepted by the same 
general requirements as in the search for secondary vertices, and not selected by the 
vertex analysis. The events were accepted as candidates if: 

• the first measured point of at least one of the tracks was within 12 cm of the beam 
spot in the plane transverse to the beam axis; 

2 The impact parameter is defined as the distance of closest approach of a charged particle to the reconstructed primary 
vertex. The impact parameters in the Rrfi and Rz plane are evaluated separately. The sign of the impact parameter is 
defined with respect to the jet direction. It is positive if the vector joining the primary vertex to the point of closest 
approach of the track is less than 90° from the direction of the jet to which the track belongs. In events with two particles, 
each reconstructed track is considered as a jet. 
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Figure 7: Acollinearity for real data (dots), a simulated signal of m fl = 60 GeV/c 2 
decaying with a mean distance of 50 cm (blank histogram), and expected simulated 
Standard Model background (cross-hatched histogram) plus cosmic background (dark 
grey histogram). This last background is normalized in order to make the first bin of 
SM background plus cosmic radiation coincide with the corresponding value of real data. 
The selection on this variable is explained in the text and is shown with an arrow. 

• both tracks were reconstructed in the TPC to guarantee a good track reconstruction 
quality; 

• at least one of the tracks had an impact parameter larger than 0.2 cm in the plane 
transverse to the beam axis, to remove standard model events; 

• the ratio of the maximum impact parameter over the minimum impact parameter 
in the Rcj) plane was smaller than -1.5 or larger than -0.5, to reject cosmic rays; 

• the acollinearity between the two tracks was larger than 10°; 

• the angle defined by the directions of the starting points of the tracks with respect 
to the the beam-spot was larger than 3°. 

The efficiencies were derived for the different f\ masses and decay lengths by applying 
the same selection to the simulated signal events. The maximum efficiency was 29.2%, 
corresponding to a mean decay length of 2.5 cm, decreasing very fast for lower decay 
lengths due to the cut on minimum impact parameter. For longer decay lengths, the ap- 
pearance of reconstructed f\ tracks in combination with the cut on the maximum amount 
of charged tracks caused the efficiency to decrease smoothly. This decrease is compen- 
sated by a rising efficiency in the search for vertices. No dependence on the f\ mass 
was found far from the kinematic limit. The losses of efficiency for f masses near the 
kinematic limit and due to initial state radiation were also considered. 

Trigger efficiencies were studied, simulating the DELPHI trigger response to the events 
selected by the vertex search and by the large impact parameter analysis, and were found 
to be around 99%. 
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No events in the real data sample were selected with the above criteria. The number 
of expected background events at y/s = 183 GeV is shown in Table |2| for the combination 
of the vertex and large impact parameter searches. 



Observed events 





Total background 


u -°°-0.12 


Z*/ 7 -^(rr)(n 7 ) 


u - u '-0.06 


Z7 7 -(ee)(n7) 


u - uy -0.03 


4-fermion (except 77) 


u - iu -0.06 


77 — > T + T~ 


u - zu -0.06 


77 — > e + e~ 


17 +u - 39 

v - 1 ' -0.05 



Table 2: The number of observed events at y/s = 183 GeV, together with the total 
number of expected background events and the expected numbers from the individual 
background sources, for both large impact parameter and secondary vertex searches com- 
bined. 



3.2.3 Small impact parameter search 

The large impact parameter search can be extended further down to mean decay 
lengths of around 0.1 cm. Charged particles were selected if their impact parameter was 
less than 10 cm in the transverse plane and less than 15 cm in the longitudinal direction 
and their polar angle between 20° and 160°. Their measured momentum was required to 
be larger than 400 MeV/c with relative error less than 100% and track length larger than 
30 cm. Any calorimetric deposit associated to a discarded charged particle was assumed 
to come from a neutral particle. 

This search was restricted to events with 2 to 4 charged particles and missing energy 
larger than 0.3 y/s. The 77 events were suppressed by requiring that the visible energy 
(E V i S ) be bigger than 0.08 y/s and the transverse missing momentum larger than 0.03 y/s. 
The polar angle of the missing momentum was required to be between 30° and 150° and 
the total energy in the forward and backward regions (E30) was required to be less than 
10% of the total visible energy in the event. 

The events were then divided into two hemispheres using the thrust axis. The total 
momentum of charged and neutral particles in each hemisphere was computed and used 
to define the events' acollinearity. Standard e + e~ — > ffi'j) processes and cosmic ray 
events were reduced by requiring the acollinearity to be greater than 10°. The charged 
particle with largest momentum in each hemisphere was selected (leading particle). The 
following quality requirements were only applied to the leading charged particles: the first 
measured point of the tracks had to be within 50 cm of the beam spot in the xy plane, 
the tracks were required to have at least a track segment beyond the ID detector, and 
away from insensitive regions of the electromagnetic calorimeter. In addition, at least 
one of the tracks was required to be reconstructed with the TPC. 

The standard e + e~ — ► //( 7 ) and cosmic backgrounds were reduced by requiring the 
angle between the leading particles in the xy plane to be less than 3 radians. 99(7) and 
4-fermion events were further rejected by requiring \Jp\ + P2 (where p\ and P2 are the 
momenta of the leading particles) to be smaller than 0.03\/s. To reduce Bhabha events 
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the total electromagnetic energy of the leading particles (E eml + E em2 ) had to be smaller 
than 0.35^/s. By requiring that any leading track with an impact parameter larger than 
1 cm in the xy plane should be reconstructed by the TPC and at least another detector, 
the residual cosmics were rejected. Finally, 7 conversion events with only two tracks were 
rejected by requiring the angle between tracks at their perigee to be greater than 1°. 

The background left after the selection described above consists mainly of events con- 
taining t pairs in the final state (7* / 'Z* — > tt and WW — > tutu). To reject these events, 
the variable b c = \fb\ + b^ was used, where b\ and 62 are the impact parameters of the 
leading particles, defined in the xy plane. Figure |8] shows the b c distribution of the se- 
lected real data, the total residual simulated background and the b c distribution of one 
set of simulated signal events, with an arbitrary scale. A cut of b c > 600/im was chosen 
to reject most of the remaining background. 




Figure 8: b c distribution for real data (dots), expected Standard Model background 
(cross-hatched histogram) and a simulated signal of m fl = 50 GeV/c 2 and a mean decay 
length of ~ 1 cm (in arbitrary scale). 

Applying these cuts to the simulated signal events, the efficiency turned out not to 
depend separately on the centre-of-mass energy and on the f\ mass but rather on the 
f\ decay length in the laboratory system, which is determined by both these variables. 
The maximum efficiency was ~ 40% for a mean decay length of ~ 2 cm. The cut on b c 
caused the efficiency to drop at small decay lengths (~ 15% at 1 mm), whilst at large 
decay lengths a loss of efficiency was due to the upper impact parameter cut used in the 
track selection. 

The non-negligible background contributions, normalized to the integrated luminosi- 
ties of the four samples, are shown in table || together with the number of selected data 
events. As expected, the main sources of background come from channels containing r's 
in the final state. 
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130 GeV + 
136 GeV 


161 GeV 


172 GeV 


183 GeV 


Observed events 










1 


Total background 


u - iJ -0.08 


o.33i8:g 


n in+u.12 

u - ±y -0.05 


iy ' -0.27 


Z*/ 7 ->(rr)(n 7 ) 


V - L J -0.08 


o.2ii8;g 


u - u '-0.04 


n qq+u.2y 

u - JJ -0.22 


77 — > T + T~ 


u - uu -o.oo 


09+ Ui4 

u - UJ -0.06 


00 +U U8 
u - uu -o.oo 


u - zu -0.12 


WW 




u - uo -0.01 


u - ±z -0.03 




ZZ 








0.02^;- 



Table 3: Expected simulated SM background events and selected data events at the 
various centre-of-mass energies for the small impact parameter search. 

4 Results and interpretation 

4.1 Neutralino pair production 

No event passes the selections at y^i = 161 or 172 GeV. Two events were observed to 
pass all the cuts at -Js = 183 GeV. One of them is shown in fig. || 

Their main features are listed in Table |j. Both of them can be interpreted as being 
4-fermion events. The first event has an electron, a pion, and two unidentified low 
momentum particles. The event could be described as 7*7*, each virtual photon going 
into a pair of rs. The second event contains a muon, two energetic electrons and a pion. 
It could be described as a Z* 7 * event, with Z* e + e~ and 7* — > t + t~. 





Candidate 1 


Candidate 2 


r 


0.14 


0.31 


Pt 


8.7 GeV/c 


9.2 GeV/c 
63.8 GeV/c 2 




139.5 GeV/c 2 


Thrust 


0.91 


0.84 


E3o/E„i S 


0.45 


0.62 


E7 1 


28.0 GeV 


78.6 GeV 


Acoplanarity 


8.6° 


15.9° 


Num. of charged particles 


4 


6 


Min. angle between jets 


63° 


26.° 


P of leading particle 


17.8 GeV/c 


43.7 GeV/c 



Table 4: Some characteristics of the two candidates found at yfs = 183 GeV. 

Since no evidence for a signal was found in the data, a limit on the production cross- 
section for neutralino pairs was derived for each (m^o,m^) combination. A statistical 
error of ±1.5% was assumed for the signal efficiency. 

In what follows, the model described in reference ||] will be used in order to derive 
limits. This is a general model which assumes only radiatively broken electroweak sym- 
metry and null trilinear couplings at the messenger scale. The corresponding parameter 
space was scanned as follows: 1 < n < 4, 5 TeV < A < 900 TeV, 1.1 < M/A < 9000, 
1.1 < tan/3 < 50, and fi > 0, where n is the number of messenger generations in the 
model, A is the ratio between the vacuum expectation values of the auxiliary component 
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Figure 9: yz view of the second of the two neutralino pair production candidates. From 
above, clockwise, the particles are identified as a muon, a pion, an electron (with associ- 
ated radiated photon) and an electron. 



superfield and the scalar component of the superfield and M is the messenger mass scale, 
tan (3 and fi are defined as for the MSSM. 

Figure 10 shows the 95% C.L. upper limit on the Xi P arr production cross-section at 
y/s = 183 GeV as a function of m^o and after combining the results of the searches 
at s/s = 161, 172 and 183 GeV with the maximum likelihood ratio method [P8"| . For 
different number of messenger generations, the ratios between production cross sections 
at different energies are bound to vary within certain limits. The same happens when 
considering scenarios with higgsino- or gaugino-like neutralinos. Figure ID presents as 
an example the case of n = 3 and gaugino-like Xi- F° r t ne other scenarios considered 
in this study (1 < n < 4, and gaugino- or higgsino-like neutralinos), the maximum 
difference with respect to figure |l^ occurs in the region where m^o < 80 GeV/c 2 and 
< 65 GeV/c 2 , and is not bigger than 10%. 

Given the aforementioned limits on the production cross-section, some sectors of the 
(m^o,m fl ) space can be excluded. In order to achieve the maximum sensitivity, the 
results from two other analyses are taken into account. The first is the search for fi pair 
production in the context of the MSSM. In the case where the MSSM Xi is massless, 
the kinematics correspond to the case of fi decaying into a r and a gravitino, except for 
spin effects, which are not taken into account in SUSYGEN. The second is the search 
for lightest neutralino pair production in the region of the mass space where Xi is the 
NLSP [Tl|] (the region above the diagonal line, i.e. > m^o). Within this zone, the 
neutralino decays into a gravitino and a photon. 

As an illustration, fig. |Tl| presents the 95% C.L. excluded areas for the case n = 2 
and gaugino-like neutralinos in the m^o vs. plane. The positive-slope dashed area is 
excluded by this analysis. The resulting 95% C.L. lower limit on the mass of the lightest 
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Figure 10: 95% C.L. upper limit of the Xi P a i r production cross-section (in picobarn) at 
\fs = 183 GeVafter combining the results of the searches at yfs = 161, 172 and 183 GeV, 
as a function of m^o and m fl for the case n = 3 and gaugino-like neutralinos, where n is 
the number of messenger generations. The diagonal and vertical lines show respectively 
the limits m^o = w,f + m T and m%o = \/s/2. 



neutralino is 78 GeV/c 2 . The negative-slope dashed area is excluded by the analysis 
searching for neutralino pair production followed by the decay Xi ~^ The point- 
hatched area is excluded by the direct search for MSSM f\ pair production ||10|| , taking 
into account the possibility of ft — Tr mixing |7| . 

For other cases, lower limits for the mass of the lightest neutralino obtained with 
this analysis are described in table [5[ In the case of n = 1 and gaugino-like lightest 
neutralino, the NLSP is always x\i an d the lower limit is derived from the search for 
acoplanar photons [|TT|. 



n 


gaugino-like Xi 


higgsino-like Xi 




(GeV/c 2 ) 


(GeV/c 2 ) 


1 


81.0 


71.0 


2 


78.0 


71.0 


3 


77.0 


49.0 


4 


78.0 


45.0 



Table 5: The 95% C.L. lower limits on m^o for eight different scenarios. When n — 1 and 
the lightest neutralino is gaugino-like, the limit comes from the search for two acoplanar 
photons. 



4.2 Stau pair production 

No candidate was observed to pass the selection of the large impact parameter and 
secondary vertex searches while the total number of background events expected was 0.63 
(0.37 on the vertex search and 0.26 on the large impact parameter search). The results 
of these analyses were combined with those of the stable heavy lepton search described 
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Figure 11: Areas excluded at 95% CL. for n — 2, gaugino-like neutralinos and tuq < 1 eV 
in the m^o vs. plane. The positive-slope dashed area is excluded by this analysis. The 

negative-slope dashed area is excluded by the search for x\ ~> jG, and the point-hatched 
area by the direct search for stau pair production in the MSSM framework. The wiggled 
curve indicating the limit from the search for neutralino pair production, is due to a rapid 
variation in the 95% CL. limit on the production cross section. 



in ||, which considers f\ decays outside the tracking devices (R > 200 cm). For very 
large f\ masses, efficiencies around 80% were obtained by the heavy lepton search. Given 
that an event could be selected both by the vertex search and by the stable heavy lepton 
search, the correlation was taken into account. 

One candidate was observed to pass the small impact parameters search, due to a 
leading track with b\ = 1.3 mm caused by an interaction with the material of the mi- 
crovertex detector. A maximum efficiency of around 40% was estimated for this search, 
and the expected SM background was 1.97 events. 

Figure [12| shows the 95% CL. upper limit on the stau pair production cross-section 
at yfs = 183 GeV after combining the results of the searches at y/s = 161, 172 and 



183 GeV with the maximum likelihood ratio method f28| . The results are presented in 
the (m^m^) plane combining the two impact parameter searches and the vertex analysis. 
The minimum upper limits achieved for a given fx were around 0.10-0.15 pb depending 
on itiq. For 13 eV/c 2 < itiq < 150 eV/c 2 and a 70.0 GeV/c 2 ¥\, a 0.15 pb limit was 
obtained. 

The upper limits on the production cross-section were used to exclude values 
as a function of ttlq combining all LEP2 energies, assuming conservatively the f\ to be 
right-handed. The vertex analysis allows the exclusion of tr masses between 70.0 and 
77.5 GeV/c 2 at 95% CL. in the range of intermediate gravitino masses (25 to 150 eV/c 2 ), 
the stable heavy lepton search covers the high gravitino mass region (over 100 eV/c 2 ), 
while the large and small impact parameter searches cover the region of low gravitino 
masses. 

Combining these results with the results of the search for MSSM f^, allows the exclu- 
sion of stau masses below 68.5 GeV/c 2 irrespective of the gravitino mass. The results are 
shown in Figure [T^, for G larger than 250 eV/c 2 (not shown in the plot) the limit was 
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Figure 12: 95% CL. upper limit of the e + e~ — > f f production cross-section (in picobarn) 
at y/s = 183 GeV after combining the results of the searches at y/s = 161, 172 and 
183 GeV. Results are shown in the (m^mfj plane. Searches for events containing 
charged tracks with small impact parameter, large impact parameter or vertices are used. 



80.0 GeV/c 2 obtained from the stable heavy lepton search |J. Following |4| as in sec- 
tion |4.1|, figure ITJ shows the 95% C.L. excluded areas for the case of n = 2, gaugino-like 



neutralinos and m 



G 



40 eV/c 2 in the m^o 



vs. mfj plane. The negative-slope dashed 



area is excluded by the analysis searching for neutralino pair production followed by the 
decay \\ ~* The point-hatched area is excluded by this search taking into account 
the possibility of Tl—Tr mixing J/J. The resulting 95% C.L. lower limit on the mass of the 
lightest neutralino is 62 GeV/c 2 , and that for the stau 60 GeV/c 2 , from the neutralino 
pair production search. 

By comparing figures 11 and 14 it can be seen that the exclusion power of isolated 
photon searches decreases as the mass of the gravitino increases. It can also be seen that 
the area excluded by the stau pair production searches increases with uiq. 



5 Summary 

Lightest neutralino and stau pair production were searched for in the context of light 
gravitino scenarios. It was assumed that the f\ is the NLSP and that the x? is the 
NNLSP. Both searches were used in order to explore the {m^m^) plane in different 
domains of the gravitino mass. 

The search for neutralino pair production produced two candidate events to be com- 
pared to 0.77±0.16 events expected from the SM background for the samples ranging 
from \fs = 161 GeV up to 183 GeV. An upper limit on the corresponding production 
cross-section between 0.2 and 0.3 pb was set at 95% C.L. in the kinematically allowed 
region. 

The search for the pair production of long lived staus produced one candidate for 
the small impact parameter method and none for the large impact parameter and vertex 
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Figure 13: Exclusion region in the (mQ,mf R ) plane at 95% CL. for the present analysis 
combined with the stable heavy lepton search and the MSSM tr search, using all LEP-2 
data. The positive-slope hatched area shows the region excluded by the small impact 
parameter search. The negative-slope hatched area shows the region excluded by the 
combination of the large impact parameter and secondary vertex searches. 



methods, whereas totals of 2. 681^3^ and 0.631^2 events were expected from the simulated 
SM background, respectively. An upper limit on the stau pair production cross-section 
was set as a function of its mass and that of the gravitino, between 0.1 and 2 pb at 95% 
CL. in the kinematically allowed region. This result, together with the search for staus 
within the MSSM framework and stable stau production, allow the DELPHI collaboration 
to set the lower limit on the mass of the tr at 68.5 GeV/c 2 at 95% C.L.. 
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